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О модели Рихарда в задачах патологических 
образований с учетом иммунного ответа 


У статті розглядається побудова математичної моделі процесу росту патологічного утворення. Розглянуто 
існуючі моделі росту клітинних популяцій та запропоновано розглядати ріст патологічних утворень 
на основі динамічної моделі Ріхарда. П розв'язок дозволяє виявити механізми, що обумовлюють 
руйнування аномальних клітин і зберігання сусідніх. 
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В статье рассматриваєтся построение математической модели процесса роста патологического образования. 
Рассмотреньт существующие модели роста клеточньжх популяций и предложено рассматривать рост 
патологических образований на основе динамической модели Рихарда. Бе решение позволяєт вьїявить 
механизмьі, вследствие которьїх разрушаются аномальнье клетки и сохраняются соседниє. 
Ключевьтге слова: патологический процесс, модель Рихарда, узловьте образования, иммунньій ответ. 


Вступ 


Останніми роками значна увага дослідників привернута до патологічних утворень 
в людському організмі, зокрема - до вузлових утворень щитоподібної залози (ЩЗ). 
Це пов'язано з кількома причинами: об'єктивним зростанням рівня захворюваності 
внаслідок радіаційних наслідків аварії на ЧАЕС, поліпшенням якості діагностики 
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вузлової патології ЩЗ завдяки широкому запровадженню сонологічного скринінгу 
та тонкоголкової аспіраційної пункційної біопсії (ТАПБ) вогнищевих уражень ЩЗ (1. 

Вузол найчастіше є пухлиною щитовидної залози (доброякісною або злоякісною). 
Досить рідко у формі вузла може бути кіста або змінена внаслідок хронічного запаль- 
ного процесу тканина щитовидної залози. Для того аби виросла пухлина, необхідно, 
щоб або протипухлинний імунітет був пригнічений, або кількість пухлинних клітин 
була дуже великою. Така ситуація може виникнути, коли щитовидна залоза зазнає 
впливу факторів, які провокують утворення пухлинних клітин (іонізуюче опромінення, 
хімічні мутогени) або стимулюють поділ клітин залози (підвищена продукція ТТГ, 
йодна недостатність). За таких умов утворюється група пухлинних клітин, яку імунна 
система не знешкодила. Ці клітини дають початок пухлині, яка росте автономно. У ро- 
ботах |2), (3| розглядають ріст пухлинних популяцій на основі динаміки Гомпертца, 
проте недостатньо вивчене питання розвитку загального патологічного утворення. 

Мета даної роботи - здійснити аналіз існуючих моделей росту патологічних 
утворень та запропонувати динамічну модель доброякісних патологічних утворень з 
урахуванням імунної відповіді. 


Аналіз моделей росту чисельності популяцій 


Останнім часом з'явилося чимало робіт, присвячених динаміці росту патоло- 
гічних утворень (2-9). Більшість запропонованих моделей мають дуже складну мате- 
матичну структуру: нелінійні диференціальні рівняння, стохастичні диференціальні 
рівняння, функціонально-диференціальні рівняння, тому їх складно досліджувати. 
Важливим залишається питання про доцільність вибору математичної моделі, яка б 
найточніше відобразила природній процес росту патологічних утворень. Отже, роз- 
глянемо найчастіше вживані математичні моделі, що описують ріст пухлини. 

Використовують такі типи моделей: Берталанфі (4, р. 89-133|, Ріхарда Г|31, логі- 
стичну (|6|, Гомпертца |2|, (3). Людвіг Берталанфі |4| запропонував загальну форму 
для моделі зміни чисельності популяції: 


ХО) ву ар (хо (1 
ай 


де х - чисельність, В,К, р,п - додатні сталі. Рівняння (1) називається узагаль- 
неною моделлю Берталанфі з параметрами р і п. Рівняння Берталанфі є узагальненням 
ряду моделей популяційної динаміки. 

Так при р-1 та п - І отримуємо логістичне рівняння (б: 


т 


ан вход 1 Я), 0) 


Іншим прикладом моделі, яка може бути отримана з моделі Берталанфі, є рів- 
няння Гомпертца |21, |ЗІ: 


ах(т) Я. 
ЕРА -сх(1)п да (3) 
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де с» - деяка стала. Рівняння (3) отримується з моделі Берталанфі при 


В - г , р- 1 та істотно малих значеннях п. А саме: 
п 
; ЕЗі (хо - 
5 ка аю а о ПР 
іа хоЗ|їо 29) і п о о НОЙ 
п-» 0 і К п-» 0 й п-» 0 і 


2 зсх(Т)Їп -ехСУПтресг)) -- ак) з еж(гУСщК| - Паре(2))) є сх(т) п 


"Р 
і 


5 
хо 
Як правило, функцією Гомпертца описується ріст пухлини |2, (ЗІ. 

Ще однією моделлю, яка випливає з моделі Берталанфі при р - 1, є модель Ріхарда 
(або узагальнена логістична) |41: 


аха) 
Фо 


п 
Вх(т) ін20) (4) 
Кк 
У роботі (2) розглянуто клас моделей на основі динаміки Гомперца, що засто- 
совувалися при дослідженні моделей пухлинного росту. Процеси росту патологічних 
утворень більш широкого класу (доброякісних або злоякісних) вимагають загальніших 
математичних моделей, до яких, як показано вище, належить модель Ріхарда. 


Постановка задачі 


Слідуючи загальній методології, запропонованій в роботі (10) для інфекційної 
імунології та розвиненій в роботах |2), |3| для моделей протипухлинного імунітету, 
будемо вважати, що основними діючими факторами процесу росту загального пато- 
логічного утворення (вузлів, фолікулярних утворень та ін.) є такі величини: 

1. Концентрація клітин патологічного утворення (1). 

2. Концентрація антитіл Е(г). Під антитілами розуміють субстрати імунної си- 
стеми, що нейтралізують рецептори клітин патологічного утворення. 

3. Концентрація плазматичних клітин С(т). Це популяція носіїв 1 виробників 
антитіл. 

4. Відносна характеристика росту патологічного утворення т(т). 

Переходимо до побудови рівнянь моделі. Перше рівняння буде описувати зміну 
чисельності клітин патологічного утворення в органі (наприклад в щитовидній залозі): 


п 


11) 


о а ЩО) іє о 
її 


-ї ЕІ). (5) 


Вираз у лівій частині рівняння (5) описує приріст клітин патологічного утворення 
4-7) за інтервал часу аї за рахунок розмноження згідно з моделлю Ріхарда (4). Цей 


приріст «Щ(т) пропорційний до є, - «народжуваності» клітин патологічного утворення. 
Член у, К(1)Щ(1) описує кількість клітин, які нейтралізують антитіла ЕК (т) за інтервал 
часу ат. Кількість таких патологічних клітин буде пропорційна як до кількості анти- 
тіл в організмі, так і до загальної кількості клітин патологічного утворення; у, - кое- 
фіцієнт, що визначає ймовірність нейтралізації (руйнування) патологічних клітини 
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антитілом. 0, - межа насичення популяції клітин патологічного утворення; п - пара- 


метр форми моделі Ріхарда, зміст якого визначається якісним характером патологіч- 
ного процесу. 

Друге рівняння описуватиме ріст плазматичних клітин. З цією метою викори- 
стаємо таку гіпотезу про формування каскадних популяцій плазматичних клітин. 
Оскільки в нашій моделі під антитілами розуміємо субстрати, що здатні зв'язуватися 
з рецепторами клітин патологічного утворення, то кількість лімфоцитів, стимульованих 
таким чином, буде пропорційна до Щ1)Е (т) . Приходимо до співвідношення, що описує 
приріст плазмоклітин над нормальним рівнем плазмоклітин в здоровому організмі: 


І ари -к)БО -л)- но (СИ)- Су), (6) 


Перший доданок у правій частині (6) описує генерацію плазмоклітин, т - час, 
протягом якого здійснюється формування каскаду плазматичних клітин,  - коефі- 
цієнт, що враховує ймовірність зустрічі антиген - антитіло, збудження каскадної 
реакції 1 кількість новоутворених клітин. Другий доданок у цій формулі описує змен- 
шення кількості плазмоклітин за рахунок старіння, й, - коефіцієнт, обернений до 


часу життя плазмоклітин. 
Для отримання третього рівняння підрахуємо баланс кількості антитіл, що реа- 
гують з рецептором клітини патологічного утворення: 


4Е 
РО си, зту, ВШ )БЕ(І). (7) 


Перший доданок справа описує генерацію антитіл плазматичними клітинами, 
другий доданок описує зменшення кількості антитіл за рахунок зв'язку з рецепторами 
клітин патологічного утворення та старіння антитіл. 

Побудова рівняння не враховує послаблення життєдіяльності організму під час 
захворювання, що пов'язане зі зменшенням активності органів, які забезпечують ви- 
роблення імунологічного матеріалу: лейкоцитів, лімфоцитів, антитіл і т.д., необхідних 
для боротьби з антигенами. Приймемо гіпотезу про те, що продуктивність таких органів 
пов'язана з розмірами пошкодження органа-мішені. З цією метою введемо в розгляд 


рівняння для відносної характеристики пошкодження органа-мішені. Нехай У -- 


п 


об'єм утворення в неушкодженому органі, а У, - відповідна характеристика патоло- 


гічного утворення. Введемо величину ПТ - ступінь ушкодження органа 
орної ов, (8) 
У 
Для утворення в неушкодженому органі - пт-0, а для органу, повністю знище- 
ного патологічним утворенням, -- т - 1. Для цієї характеристики пропонується рівняння: 
ВЕ РИ (9) 
аї 
У рівнянні (9) перший доданок описує необоротні зміни в клітинах патологічного 
утворення, а другий - відновлення утворення за рахунок регенерації тканини. При 
сильному ушкодженні життєво важливих органів продуктивність вироблення анти- 
тіл падає. В нашій моделі фактор ушкодження життєво важливих органів можна вра- 
хувати в рівнянні (6), замінивши коефіцієнт а на добуток а (т), де 
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де К - стала, що залежить від локалізації захворювання. В інтервалі 0 т т! 
функція 25 (т) дорівнює одиниці. Це означає, що функціонування імунологічних органів 


в цьому інтервалі не залежить від складності захворювання. Але при т «т 1 їх 
продуктивність швидко падає. 

Отже, підсумовуючи вищенаведені викладки, приходимо до таких біологічних 
припущень: 

1. Популяції клітин патологічного утворення, антитіл, плазматичних клітин є 
однорідними. 

2. Зміна чисельності популяції клітин патологічного утворення підлягає законам 
динаміки Ріхарда. 

3. Їмунна реакція індукується рецепторами клітин патологічного утворення та 
антитілами одного виду. 

4. Концентрація рецепторів в момент часу г пропорційна до кількості клітин 
патологічного утворення ЩІ). 


5. Клітини патологічного утворення мають стримуючий вплив на зростання 
популяції антитіл. 

Таким чином при відсутності жодного лікувального впливу приходимо до такої 
системи диференціальних рівнянь: 


п Де що -у БОЦІ), 10) 
а 8, 
пе -5(таї от) -т)- мо (С()-С,), (1 
- 9 -Б,СІ)- (и, ту ШОУЕОІ), (12) 
о ОЛИ СЦІ)- у, т(г), о 
7 


Тут ЩІ) - кількість клітин патологічного утворення, С(т) - концентрація плаз- 
матичних клітин, Е(7) - концентрація специфічних антитіл, т(т) - ступінь ушко- 
дження органа, у, - коефіцієнт, що визначає ймовірність нейтралізації (руйнування) 
клітини патологічного утворення антитілом, Р, - швидкість виробництва антитіл 
однією плазматичною клітиною, а, - коефіцієнт, що зумовлює ймовірність зустрічі 
антиген-антитіло, и, - коефіцієнт, обернений до часу життя плазматичних клітин, 
и, - коефіцієнт, обернено пропорційний до часу розпаду антитіл, у, - коефіцієнт, 
що враховує швидкість відновлення пошкодженого органу, п - число специфічних 
антитіл, що потрібно для нейтралізації одного антигена, с - коефіцієнт, що визначає 
швидкість загибелі клітин за рахунок пошкоджуючої дії антигенів, 5 (т) - неперервна 
незростаюча функція (0 «4 2 (т) 4 1), що характеризує порушення нормального функ- 
ціонування імунної системи через значне пошкодження органа-мішені. 

Перераховані параметри є додатні 1 специфічні як для виду рецептора, так і для 
органа 1 для конкретного організму. 

При цьому ЩІ), С(1), Е(1), т(Т) є С! рав оо). 

Задано неперервні початкові умови на 1 є Їх, еру 

щі) з 1), СИ) З СО), піт) сть). а4 
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Зазначимо, що система типу (10) - (13) є узагальненням відомої моделі Г.І. Мар- 


чука, де рівняння (10) - логістичного типу. 


Приклад 


За допомогою комп'ютерної програми «Інтегратор систем з післядією (ІХТЕСВА- 


РОЗТ)» (Свідоцтво про реєстрацію авторського права на твір ХМ» 6027, 30.07.2002, 
Державний департамент інтелектуальної власності МОН України) проведено чисельний 
аналіз моделі (11) - (15). Розглянуті такі значення параметрів: 


а, 2 0,00396, 0, -14, у, - 0,008, с - 107, п, 20,5, ре017, ш, 2017, п -10, 


и, 7012, 5. «0,7 «175, п-8, 


10 «т «0,1 
е(т) з 
а-ту/(а0/9)01«тх«ї1 


При т є|-т,0) справедливі такі початкові умови: Щі) - 1, С(1)- 1, К(7)- 1, 


т(т) - 0. На рис. І показано розв'язки рівнянь моделі. 


1.75 9.7130.21.176 32.639 44.102 55.565 67.028 78.49189.954 101.41 112.88 124.34 135.80 147.26 158.73 170.19 181.65 193.12 


Рисунок 1 - Чисельне моделювання системи (11) - (15). На графіках показано 
1-1000Г(), 2-С()/100, 3-К()/80, 4-т(т) 


Метою подальших досліджень повинно стати обчислення станів рівноваги та 


отримання умов асимптотичної стійкості в класах рівнянь популяційної динаміки. 


Література 


1. 


Шідловський В.О. Ефективність лазеріндукованої інтерстиціальної термотерапії у лікуванні вузло- 
вого зоба / В.О. Шідловський, О.В. Шідловський, О.І. Карел // Український журнал хірургії. - 
2009. -- Хо 4.- С. 137-139. 

Марценюк В.П. Построение и изучение устойчивости модели противоопухолевого иммунитета / 
В.П. Марценюк // Кибернетика и системнь(й анализ. - 2004. - Моб. - С. 123-130. 

Марценюк В.П. Об устойчивости в модели иммунной защитьт с учетом нарушения функциониро- 
вания органа-мишени: метод вьтрожденньх функционалов Ляпунова / В.П. Марценюк // Киберне- 
тика и системньтй анализ. - 2004. -- Мо 1. 

Ва|гег 7. Мафетайса! подейпе ої (штог 5гом'Ю Кіпейсз / Ва)гег 7., МиК-Раміоміс 5., Нигак М. - 
Возіоп, 1997. 


272 «Искусственньюй интеллект» 172012 


Про модель Ріхарда в задачах росту патологічних утворень... вм 


5. Мафетайсаї подеШпє ої шштог ягомір / Р. Стеси, А.5. Сагягеа, А.Т. Стеси, А. Мі5еїпезси // Котапіап 
Керогіз їп РПузіс5. - 2007. - Мої. 59, Хе 2, - Р. 447-455. 

6. М/ізе 5.М. Тргее-дітепзіопа! пиіцйзресіе5 попійпеаг (штог гор. Модеї! апа питегісаї плефоай / 5.М. М/іке, 
7.5. помепетир // 7. ТПеогеї. Віої. - 2008. - Мої. 253. -Р. 524-543. 

7. Сегі5сп А. Мафетайса! тодеШпе ої сапсег сеї іпуазіоп ої пязие: Г.осаї апа поп-Іоса! тоаеія апа Ше 
еНесі ої адрезіоп / А. Сегізсі, А... Сраріаїп // 7. ТПеогеї. Віої. - 2008. - Мо 250. - Р. 684-704. 

8. Адат ).А. А 8игуеу ої подеїз Їог (штог-іпатипе 5у5іет дупатіся / 7.А. Адат, ХХ. ВеПото, Вігкрайеег. - 
Возіоп, 1997. - С. 389-133. 

9. Мафетайса! птодейпе ої сапсег гайіоуігоШегару / Р. Діпеї, М.Ю. Сазсіпо, К. )о5і.с, 5.). Ки5заеії, 7. 
Ва)угег // Маю Віозсі. - 2006. - Мої. 19.-Р. 55-78. 

10. Марчук Г.Й. Математические модели в иммунологии / Марчук Г.И. - М. : Наука, Главная редакция 
физико-математической литературьі, 1980. - 264 с. 


І леегагига 


1. 5рідоузє'Кі) М.О. ОкКгаїпз5'Кі) г2бигпа! рігогеїї. 2009. Хе 4. 5. 137-139. 

2. МагсепуиКк У.Р. Кібегпейка 1 515їепапу| апаП7. 2004. Мо 5. 5. 123-130. 

3. МагсепіиК У.Р. Кібегпе-йКа і 515їепапу) апа. 2004. Хо 1. 

4. Ва)гег 7. Мафетайса! тодеПпе ої іштог єгом'Б Кіпейсз. Возіоп. 1997. 

5. Стеси Б. А.5. Котапіап Керогіє їп РПузісз. 2007. Мої, 59, Хе 2. Р. 447-455. 

6. Мізе 5.М. 7. Треогеї. Віої. 2008. Мої. 253. Р. 524-543. 

7. Сегіз8сі А. ). Треогеї. Віої. 2008. Хо 250. Р. 684-704. 

8. Адат 7.А. А 5игуеу ої тодеї8 Гог (итог-іпатипе 5убіет дупатіся. Возіоп. 1997. 5. 89-133. 

9. Діпєй р. Май Віовсі. 2006. Мої. 19. Р. 55-78. 

10. МагепикК С.І. Магетайсрезкіе подає у ітпипоїогії. М.: Кашйка. СіІаупа|а гедаКксіа ПліКо- 
тагтайспезко) Піегаїшгу. 1980. 264 5. 


У,Р. Магі5епуик, О.А. Ва»гі)-Хауаїя 
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ОЇ те Іттипе Кезроп5е 

Ехізцпє піодеіє ої сей рориіайоп ягомії аге апаїугед. Боиг еппрігіса! піодеї5 
демеїоредй Гог Біоіоєісаї 5укіетя ягом'ї, 5исП аз Вегіаїапійу, Кіспага, Іоєб15йїс, Сопрегі7 
аге сопзідегей. Тре ацПог5 рге5епіед декегтіпі5йс пзаїретайса! подеї8, муБісП аге пло5і 
ойеп и5ей Їог (штог ягоу/П дебсгірпоп. ІП 15 ргоро5ед іо соп5ідег єгом/Ю ої раїоїіогісаї 
Гогтабоп оп Ба5і5 ої Кіспагда5 дДупатіс подеі. 

Такіпє їпіо ассойпі Фе іаптипе гезропзе плоаеі, едиайопя ої а плодеі аге сопзігисіва. 
Тре Нт5і едиайоп дезсгібез Ше спапее5 ої се!! питбег ої рафоіодіса! Їогтайоп їп а питап 
роду (Гог ехапріе, їп Шугоїд). ТПе 5есопа едиайоп декстібез ріазта сеїя стгому р. Тре Фіта 
едиайоп дезсгібез Ше срапее5 ої питбег ої апійбодіе5, у/бісП геасі Бу (Фе гесеріог ої 
раФЧоїоєвіса! се!!8. ТПе Гоц'П едиайоп дебсгібез Ше ехіепі ої огсап датаєе. 

Тре пафетайсаї птоадеі Їог єгом' ргосе58 ої рафоіоєісаї! Гогтайоп (аКіпе, іпіо 
ассоцпі Ше іттипе гезроп5е 15 биї. Те пиплегіса! ехатріє ої а плодеї мії 50те 
ріоіовісаї аззитріопея 15 ргезепіед. 


Стаття надійшла до редакції 11.11.2012. 
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